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We explore in depth the validity of a recently proposed scaling law for earthquake inter- 
event time distributions in the case of the Southern California, using the waveform 
cross-correlation catalog of Shearer et al. Two statistical tests are used: on the one 
hand, the standard two-sample Kolmogorov-Smirnov test is in agreement with the scaling 
of the distributions. On the other hand, the one-sample Kolmogorov-Smirnov statistic 
complemented with Monte Carlo simulation of the inter-event times, as done by Clauset 
et al., supports the validity of the gamma distribution as a simple model of the scaling 
function appearing on the scaling law, for rescaled inter-event times above 0.01, except 
for the largest data set (magnitude greater than 2). A discussion of these results is 
provided. 
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1. Introduction 

In the last years considerable attention has been addressed to the distribution of 
inter-event times in natur al hazard s, in particular ear 

thquakejl| 2|3|4|5|617|8|9110|^ 

but also in human response J^^ESIISlg^j^jj^ social 

behavioil ^4 | 15n(j | 17 : 

In some of these 

systems, the shape of the inter-event-time distribution for events above a certain 
threshold in size is independent on the threshold. Indeed, let r denote the inter- 
event time, defined as the time between consecutive events above a size threshold s 
and let Ds{t) be its probability density, then we can write 

DAT)^Rsf{Rsr), (1) 

where / is a scaling function that provides the shape of Ds (t) and Rs is the occur- 
rence rate for events above s, providing the scale of Ds{t) (and Rj^ provides the 
scale for r). 

If the size distribution follows a power law (which is not always the cas^^^, 
then Rs cx (where /3 is the exponent of the cumulative size distribution), and 
then 

i^,(r)=/(r/.^)//, (2) 
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which turns out to be a scaUng law, equivalent to those obtained in the study 
of critical phenomenal!^. The law reflects a scale-invariant condition: there exist 
a change of scale in r and s (a linear transformation) that does not lead to any 
change in the statistical properties of the process, at least regarding the inter- 
event-time probability density. The function / is just the scaling function /, except 
for proportionality constants. 

Why is this scaling law of some relevance or interest? In general, when events are 
removed from a point process (as it is done in our case by raising the size threshold), 
the resulting inter-event-time distribution changes with respect to the original one, 
and a scaling law as Eq. H]) does not apply. However, for high enough thresholds s 
(for the extreme events that are of interest in hazard assessment studies), and when 
the events are randomly removed, it is expected that the resulting time process 
tends to a Poisson process (which means that the occurrence of the extreme events 
is independent on the history of the process, and just some multi-faced dice, thrown 
in continuous time, decides if an event takes place or not). From the point of view 
of statistical physics, the Poisson process constitutes a trivial fixed-point solution of 
the renormalization equations describing the thi nning o r decimation performed in 
event occurrence when the size threshold is raisecP^'^QI Yot event occurrence on a 
large spatial scale, as for instance worldwide earthquakes, there is a second reason to 
expect exponential inter-event-time distributions: the pooled output of several time 
processes (i.e., China seismicity, superimposed to Jap an seismicity, etc..) tends to 
a Poisson process if the processes are independent!^. 

It is therefore surprising not only that the scaling function / is not exponential, 
but also that a non-exponential scaling function exists. In the case of earthquakes^ 
(and fractures^ '^) / is approximated by the so called gamma distribution, with 
parameters 7 and a, 

fix) = V (-)'"%-^/°, for 2: = i?,T > m, (3) 
aT{'j,m/a) \x/ 

where r(7, m/a) is the complement of the incomplete gamma function (not normal- 
ized), r(7, u) = u^^^e^^''du. The cutoff value m is not considered a free param- 
eter but fixed and the scale parameter a is not independent but can be obtained 
from the value of 7 and m taking into account that [x) = [Rst) = xf{x)dx = 1 
(using that Rg is the inverse of the mean inter-event time). For stationary seismicity, 
as well as from fracture and nanofracture experiments, the shape parameter 7 turns 
out to be close to 0.7, see Refs. IH IHl [3 The reason to disregard a;— values below m 
is due, on the one hand, to the incompleteness of seismic catalogs on the shortest 
time scales and to the existence errors in the determination of the inter-event times 
when these are small, and on the other hand to the breakdown of the stationarity 
condition in those short time scales by small aftershock sequences. 

The usual way to establish the validity of a scaling law such as Eq. ([T]) is by 
plotting the different rescaled quantities together (in our case inter-event-time dis- 
tributions for different thresholds) and judge visually if they collapse onto a single 
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curve or not. It would be nice if one could put some numbers into the quality of 
the scaling and the fit of the scaling function and test their limits of validity. Let us 
note that Kagan has argued that one of the reasons because theoretical physics has 
failed not only to predict but to explain earthquake occurrence is due to the poor 
use of statistics by the researchers in the field^. Indeed, "the quality of current 
earthquake data statistical analysis is low. Since little or no study of random and 
systematic errors is performed, most published statistical results are artifacts." We 
believe this criticism has applicability beyond the case of statistical seismology. 

In this paper we will first use the Kolmogorov-Smirnov two-sample test in or- 
der to evaluate the fulfillment of the inter-event-time scaling law ^ in Southern 
California seismicity. Next, the goodness of the fit of the scaling function ([3]) to 
the rescaled inter-event-time densities will be tested by adapting the procedure 
introduced by Clauset et al. , consisting in maximum likelihood estimation of pa- 
rameters, Kolmogorov-Smirnov one-sample statistic evaluation, and Monte Carlo 
simulation of the inter-event times in order to compute the distribution of the 
statistic. 



2. Data 

The seismological data used will be the Southern California waveform cross- 
correlation catalog of Shearer et aL^^ (for which, as far as the author knows, no 
study has published plain inter-event-time distributions Ds{t), nevertheless, see 
also Refs. [531 US)- The catalog spans the years 1984-2002 (included), containing 
77034 earthquakes with magnitude M >2. Notice that we will use M as a measure 
of size, although by the Gutenberg-Richter law it is not power-law distributed but 
exponentially distributed^. In order to recover a power-law distribution one has 
to deal with the seismic moment, or the energy, which are exponential functions of 
the magnitude. 

We will concentrate in earthquake occurrence under stationary conditions. It is 
well know that earthquakes trigger more earthquakes with a rate that changes in 
time following the Omori law^H. In general, this breaks stationarity, as the rate of 
occurrence is not constant in time; however, at (relatively) large scales the resulting 
superposition of time- varying rates yields a constant rate, as it happens in worldwide 
seismic occurrence, and also for Southern California in certain time periods in which 
the largest earthquakes do not occur, see Fig. 1 of Ref.[28l Precisely for this reason, 
inter-event-times for stationary seismicity are more reliable than for non-stationary 
periods, as the l arge earthquakes present in the latter case prevent the detection 
of the small onea^, which has dramatic consequences in the computation of the 
inter-event times. 

The stationary time periods under consideration in this paper are (refining those 
in Ref. EH following Ref. 1984 - 1986.5, 1990.3 - 1992.1, 1994.6 - 1995.6, 
1996.1-1996.5, 1997-1997.6, 1997.75-1998.15, 1998.25-1999.35, 2000.55-2000.8, 
2000.9 - 2001.25, 2001.6 - 2002, 2002.5 - 2003, where time is measured in years, 1 
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year = 365.25 days (every 4 years an integer value in years corresponds to the true 
starting of the year) . 

3. Testing Scaling 

A simple way to quantify the validity of the scaling hypothesis in probability distri- 
butions can be obtained from the two-sample Kolmogorov-Smirnov (KS) test, which 
compares two empirical distributions. The procedure begins with the calculation of 
the maximum difference, in absolute value, between the rescaled cumulative distri- 
butions of the two data sets, i.e., 



as we have more than two data sets, we label them with indices k and I. The empir- 
ical cumulative distribution functions Pk (x) are calculated as the fraction of obser- 
vations in data set k below value x, and constitute an estimation of the theoretical 
cumulative distribution function, Fk{x) = Prob[variable < x] — J^Dk{x)dx. 

Obviously, the difference dki is randomly distributed, and therefore we can refer 
to it as a statistic. The key element of the KS test is that when the data sets k and 
I come indeed from the same underlying distribution F{x) = Fk{x) = Fi{x), the 
distribution of the KS statistic du turns out to be independent on the form of F{x) 
and can be easily computed. Therefore, the resulting value of dui can be considered 
as small or large by comparison with its theoretical distribution. Under the null 
hypothesis that both data sets come from the same distribution, the probability that 
the KS statistic is larger than the obtained empirical value dki gives the so called 
p— value, which constitutes the probability of making an error if the null hypothesis 
is rejected. The formulas fo r th e probability distribution of dki are simple enough 
and are given by Press et a/.'^, depending only on the number of data N in each 
of the sets; so, approximately, for large iVe, 



with Q a decreasing function taking values between 1 and (see Ref. I5T|| and 
Ne an effective number of data points (the "reduced" number, or one half of the 
harmonic mean of the number of data). Nevertheless, in order to calculate the 

value it is simpler to use the numerical routines provided in the same referencJ^ 
(in particular, the routine called probks). 

Notice that we have to compare the distributions of seismicity for M > Mk and 
M > Ml after rescaling, i.e., as a function of RkT and i?;r, respectively (otherwise, 
without rescaling, the distributions cannot be the same). In order to do that, for 
each data set, we first calculate the mean value of the inter-event time, (t)^ = -R^^, 
and then, we disregard inter-event time values such that RkT < m. The elimination 
of the smallest values increases the mean value of the remaining rescaled inter- 
event times, so, we repeat the procedure: we recalculate the mean inter-event time 
and rescale again the data by the new rate, disregarding those values below m. 



dki = rRax\/x\Pk{x) - Pi{x)\; 



(4) 



p = Prob [ KS statistic > d]= Q{WNe + 0.12 + 0.11/ ^yNe]d), 



(5) 
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The resulting data set has a mean value very close to one. We will assume that 
this procedure does not invalidate the applicability of the formulas we use for the 
calculation of the p— value. 

The rescaled inter-event-time cumulative distributions for different magnitude 
ranges are shown at Fig. [U ranging from M > 2 to M > 4, fixing m = 0.01. 
The scaling seems rather good, except for the case M > 4. Table [1] shows the KS 
statistic for each pair of distributions, as well as the corresponding p— values. Due to 
their high values (in all cases larger than 0.18 but in some others larger than 0.95), 
the null hypothesis cannot be rejected and each pair of data sets are compatible 
with the same underlying distribution, and therefore we have to agree with the 
scaling hypothesis (within statistical significance). Let us note that the p— value, 
being itself originated by a random set, is a random quantity (when different data 
sets are considered), and it turns out that the distribution of p is uniform, between 
and 1. So, there is no reason to prefer p ~ 0.9 in front oi p = 0.2. Only small 
enough values of p should lead to the rejection of the null hypothesis. 

The results for the same data using m = 0.001 (which increases A^), also shown 
in Table [11 are again in concordance with the scaling hypothesis, being the smallest 
p— value for this case larger than 0.14. Even for m = 10~^ all the p— values are 
above 0.2, except for some of the pairs involving the set with M > 2. The behavior 
of du when m is changed, which, following Ref. 23 should be a guide to chose the 
value of m (although in a different type of test, see next section), is not clear in this 
case. 

4. Testing the Scaling Function 

A different statistical test regards the goodness of the fit applied to some data. 
For instance, we can ask whether Eq. ([3]) is a good approximation to the empirical 
scaled distributions of inter-event times. Here, we will adapt the method of Ref. |23] 
to the kind of distributions that we are interested in. 

First, a fit has to be performed. A usual way of proceed in the case of long-tailed 
distributions is to minimize the squared differences between the empirical density 
and the theoretical density in logarithmic scale; however, this method shows some 
problems and involves the arbitrary estimation of the density; other problems arises 
if one fits the cumulative distributiorPSI. In contrast, maximum likelihood estimation 
avoids these difficulties by working directly with the "raw" data. 

In order to be more general, let us consider the distribution given by the prob- 
ability density. 



which constitutes the so called generalized gamma distribution, with shape param- 
eters 7 and b and scale parameter a. We consider 7 and b greater than zero, the 
opposite case can be considered as well but then the function F has to be replaced 
by its complementary function (and multiplied by -1, as (5 < 0). The cutoff value m 
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could be fixed to zero, but, as we have mentioned, for our data it is convenient to 
consider m > 0. 

The n— th moment of the distribution is given by 



7 

6 ' "o^ 



(7) 



for 7 > and S > 0. Notice that a particular case is given by the scaling function 
f{x) appearing in Eq. ([1]), for which (x) = 1, and only two of the three parameters 
are free; nevertheless we will not make use of that restriction for estimating the 
parameters. 

The method of maximum likelihood estimation is based on the calculation of the 
likelihood function L, see Ref. (23l This is given by (or, in order to avoid dimensional 
problems, proportional to) the probability per unit of x'^ that the data set comes 
from a particular distribution, given the values of its parameters i.e., 

ri X \ Prob[xi,X2, ■ ..xn^^S, a] TT \ x \ fQ\ 



where N is the number of data and we make explicit the dependence of the prob- 
ability density on its parameters. The last step assumes that each value Xi is inde- 
pendent on the rest. Naturally, this is not always the case (we know that earthquake 

09 00 00 

inter-event times are correlated ) and then the maximum likelihood method 

provides an estimation of the distribution that generates the dataset in consider- 
ation but it may be that the dataset is not representative of the process we are 
studying (due to correlations, the phase space may not be evenly sampled). 

It is more practical to work with the log-likelihood, £, which is the logarithm of 
the Hkelihood; dividing also by N, 

ii^,6,a) ^ ^IlI^hA^ = lf]lni?(..|7,<5,a), (9) 

1=1 

which, notice, can be understood as a kind of estimator of the entropy of the distribu- 
tion from the available data (with a missing —1 sign). In the case of the generalized 
gamma distribution ([6]) it is easy to get that 

>„.->„r(i(^)').„„e_(M)', a„, 

where we have omitted a term — In G that is independent on the parameters of the 
distribution, and we have introduced G as the geometric mean of the data. In G = 
{J2 lna:i)/iV, and A{S) as what we may call the ^— power mean, A{S) = {/^^'xfjN 
(which, in contrast to G, depends on the value of the parameter S] for instance, for 
S = 1, A is the arithmetic mean, but for (5 = —1, ^ is the harmonic mean). 

The best estimate of the parameters would be that that maximizes the likelihood, 
or, equivalently, the log-likelihood. The previous expression is too complicated to 
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be maximized analytically, and it is too complicated to differentiate even (in fact, 
we would need to compute the derivative of the incomplete gamma function). So, 
we will perform a direct numerical maximization (in particular we will use the 
numerical routine amoeba from Ref. '31', the function F can be computed from the 
same source using routines gammq and gcunmln) . 

Fixing (5 = 1, from which we recover Eq. ([3]) as a model of the distribution (which 
yields only one free parameter and has the advantage of being compatible with a 
Poisson process in the limit of long times), the resulting values of the parameters 7 
and a, obtained from maximum likelihood estimation, are given in Table [21 In all 
cases, except for M > 4, and if the cutoff m is not too small, the values of 7 are 
close to 0.7. 

Once we have obtained the estimators of the parameters, we can ask about their 
meaning. Maximum-likelihood estimation does not mean that it is likely that the 
data comes from the proposed theoretical distribution, with those parameters. In 
fact, maximum likelihood can be minimum unlikelihood, i.e., we are taking the less 
bad option among those provided by the a priori assumed probability model. In 
order to address this issue it is necessary to perform a goodness-of-fit test. 

Following Ref. 23: we can employ again the Kolmogorov-Smirnov test, this time 
for one sample. The KS statistic is, similarly as before 



where P{x) is the empirical cumulative distribution of the data, defined in the 
previous section, and F{x) is the theoretical proposal. For the distribution of Eq. 



The resulting values of d for our problem are also shown in Table [2l Now we can 
apply the recipe of Clauset et al. in order to select the most appropriate value of 
the cutoff m, which consists in selecting the value which minimizes d. Comparing 
between 0.003, 0.01, and 0.03, it seems clear that we should chose m = 0.01. 

At this point we could proceed as in the previous section, using the formulas 
for the distribution of d. However, that only would be right if we were not esti- 
mating F{x) from the data (if we were comparing with a theory free of parameters 
for instance). In order to know the distribution of the statistic d when the data 
are generated by the model with the parameters obtained by maximum likelihood 
estimation, we will use Monte Carlo simulations. Indeed, generating data from the 
theoretical distribution, we can repeat the whole process to obtain the statistical be- 
havior of d when the null hypothesis is true (when the data come from the proposed 
theoretical distribution), and we can do it many times, in order to get significant 
statistics. 

Schematically, the process for the calculation of p consists of the multiple itera- 
tion of the following steps: 



d = uYa,yi\jx\P{x) — F{x)\ 



(11) 



m 




(12) 
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(1) Simulate synthetic data s from the distribution given by Eq. ([3]) using the 
parameters 7 and a obtained before for the empirical data 

(2) Estimate the parameters 7^ and by fitting the synthetic data s to Eq. ([3]) 
(proceeding in the same way as described above for the empirical data, see Eq. 
(fTO]) and so on). 

(3) Evaluate the KS statistic for the distribution of synthetic data s [generated in 
(1) with parameters 7 and a] and the theoretical distribution with parameters 
7s and Ug [calculated in (2)], i.e., = maxvx|-Ps(a;|7, o) — F(x\^s,cis)\ 



We will obtain synthetic inter-event times from the gamma distribution by gen- 
erating a table of the cumulative distribution. As the probability of an event has to 
be the same independently on the random variable we assign to that event, then, 
u = F(x), where m is a uniform random number between zero and one and is also 
its own cumulative distribution. We can calculate numerically the function F{x) 
(thanks to the numerical recipes gammq and gammliiSIl) ^ but are unable to calculate 
its inverse (at least, in a reasonable computer time), so we will tabulate the values 
of F{x), for selected values of x in log scale (this is to deal with the multiple time 
scales that appear in the process, described by Eq. ([3]) or (fT^ when 7 < 1 and 
m ^ 1). To be concrete, x{k) = me"'', where fc = 0, 1, . . . and a is just a constant. 
Then, when a uniform value u is generated we can obtain the corresponding value 
of X by looking at the table and interpolating (or extrapolating) using the closest 
values of u{k) = F{x{k)). 

For the case of our interest, the p— values calculated in this way, using 1000 
randomly generated samples (which yield an uncertainty of about 3 % in p), are 
included in Tabled Taking m — 0.01 (the value arising by the application of Clauset 
et al. recommendation), we cannot reject the hypothesis that the data set comes 
from the theoretical distribution with maximum likelihood parameters, except for 
M > 2, which yields p = 0.032, which is beyond the usual onset of acceptance 
of the null hypothesis, p = 0.05. Figure [2] illustrates the reason of the rejection. 
Indeed, although the theoretical distribution is very close to the empirical one, the 
difference is large enough for the high number of data involved. Although Eq. ([5]) 
is not valid in this case, we can use it as an approximation and see how for large 
N the statistic d scales as 1/\/N {Ne — N here). As the mode of the distribution 
Q in Eq. ^ is around 0.735 and practically all the probability is contained below 
2, this means that we can expect d < 2/^/N. So, for large N, d tends to zero, and 
the KS test is able to detect any small difference between the proposed theoretical 
distribution and the "true" distribution. This means that the test is not adequate if 
we are just interested to find an approximation to the true distribution, as only the 
"true" distribution is not rejected for a sufficient number of data. For comparison, 
we show in Table [3] the results for an exponential scaling function, which is clearly 
rejected except for M > 4. 
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5. Discussion 

As another alternative, note that we have tested separately the validity of the 

scaling law and the adequacy of the scaling function given by Eq. We could 

take advantage of the scaling behavior to fit and test the goodness of fit of the 

scaling function. For instance, we could combine all rescaled data sets (for all values 

of the minimum magnitude) and proceed as in Sect. 3 for this combined data set. 

The problem is that, by virtue of the Gutenberg-Richter law, when the minimum 

magnitude is raised in one unit, the number of events decreases by a factor 10, and 

therefore, data sets with large minimum magnitudes are under-represented. Perhaps 

we could just truncate the samples in order that all of them had the same number 

of data, but that would lead to a tremendous wasting of information. 

The surprising character of the scaling law ([1]) when the scaling function is riot 

exponential has lead to some criticisms by MolcharP^I and Saichev and SornettJ^Hl. 

The latter authors propose that, for the so called ETAS model, the scaling law is 

not valid, and one has a very slow variation of the inter-event-time distribution 

when the magnitude threshold is raised. Although we have tested that the scaling 

law is consistent with the data within statistical significance, this does not mean 

that we should reject Saichev and Sornette result. Nevertheless, the simplicity of the 

scaling hypothesis makes it the most adequate model for seismicity, at least as a null 

model to contrast with other hypothesis. On the other hand, other seismicity models 

have been recently proposed, which, in contrast to the ETAS model, are fully scale 

invariant and one would expect that are characterized by scaling inter-event-time 
distributions^e 37 38 39 40^ 

Saichev and Sornette als o proy ide a pseudo-scaling function to which inter-event- 
time distributions can be fii}34i35|. 

In principle, the very same procedure used in our 
paper can be applied directly in order to fit the parameters of the Saichev-Sornette 
function and test the goodness-of-fit of the outcome. It is expected that the use of 
this new function, which has more parameters than Eq. ([3]) and models better the 
left tail of the distribution, could lead to the reduction of the cutoff m above which 
the functions are fit. We leave this task for future research. 
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Fig. 2. (a) Rescaled inter-event-time cumulative distribution for Southern California stationary 
seismicity with M > 2, fixing minimum x— value m = 0.01, together with the fit obtained by 
maximum likelihood estimation. The p— value corresponding to the KS statistic, determined by 
Monte Carlo simulation turns out to be as small as 0.032, although the fit is visually acceptable, 
(b) Same as before, in log-log scale, together with a pure power law with the same exponent, (c) 
Difference between the distribution and its fit, P{x) — Fix), which yields the KS statistic when 
its absolute value is maximized, (d) The corresponding probability density, estimated with 5 bins 
per decade, for comparison. Also the best fit and a pure power law are shown. 
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Table 1. KS statistic d (below the diagonal) and corresponding p— value 
(above diagonal, in percentage) for rcscalcd Southern- California stationary— 
seismicity inter-event-time distributions, with lower cutoffs m = 0.01 (top) 
and m = 0.001 (bottom). The scaling hypothesis cannot be rejected. 



m = 0.01 


N 


M > 2.0 


M > 2.5 


M > 3.0 


M > 3.5 


M > 4.0 


M > 2.0 


18870 




26.1% 


36.0% 


97.3% 


21.3% 


M > 2.5 


4953 


0.016 




63.2% 


86.3% 


18.3% 


M > 3.0 


1184 


0.028 


0.024 




95.7% 


22.8% 


M > 3.5 


309 


0.028 


0.035 


0.032 




20.8% 


I\l 1.0 


70 


0.125 


0.129 


n.i2G 


0.13S 






m = 0.001 


N 


M > 2.0 


M > 2.5 


M > 3.0 


M > 3.5 


M > 4.0 


M > 2.0 


19821 




41.6% 


43.8% 


14.3% 


31.9% 


M > 2.5 


5187 


0.014 




57.1% 


29.6% 


31.1% 


M > 3.0 


1268 


0.025 


0.024 




68.0% 


33.5% 


M > 3.5 


340 


0.062 


0.054 


0.044 




28.1% 


M > 4.0 


76 


0.108 


0.110 


0.110 


0.124 





Table 2. Maximum likelihood parameters 7 and a, 
KS statistic d and corresponding p— value (in per- 
centage, determined by Monte Carlo simulation) for 

rescaled Southern-California stationary-scismicity in- 
ter-event-time distributions, using several values of the 
minimum value m. 



m = 0.03 


N 


7 


a 


d 


P 


—value 


M > 2.0 


18009 


0.68 


1.35 


0.008 




1.2% 


M > 2.5 


4669 


0.67 


1.38 


0.007 




84.0% 


M > 3.0 


1122 


0.73 


1.29 


0.021 




25.7% 


M > 3.5 


287 


0.79 


1.22 


0.034 




55.7% 


M > 4.0 


69 


0.89 


1.07 


0.089 




16.3% 






m = 0.01 


N 


7 


a 


d 


P 


—value 


M > 2.0 


18870 


0.68 


1.41 


0.007 




3.2% 


M > 2.5 


4953 


0.64 


1.50 


0.009 




37.3% 


M > 3.0 


1184 


0.69 


1.41 


0.021 




21.8% 


M > 3.5 


309 


0.67 


1.45 


0.029 




73.6% 


M > 4.0 


70 


0.95 


1.05 


0.082 




29.3% 






m = 0.003 


N 


7 


a 


d 


P 


-value 


M > 2.0 


19466 


0.65 


1.51 


0.009 




0.0% 


M > 2.5 


5102 


0.62 


1.57 


0.011 




11.2% 


M > 3.0 


1247 


0.59 


1.65 


0.029 




1.4% 


M > 3.5 


328 


0.56 


1.74 


0.046 




9.0% 


M > 4.0 


74 


0.72 


1.37 


0.103 




5.8% 
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Table 3. Same as the previous table, for the ex- 
ponential distribution (7 s 1). 



m = 0.01 


N 


7 


a 


d 


p— value 


M > 2.0 


18870 


1 


0.99 


0.072 


0.0% 


M > 2.5 


4953 


1 


0.99 


0.084 


0.0% 


M > 3.0 


1184 


1 


1.00 


0.077 


0.0% 


M > 3.5 


309 


1 


1.00 


0.079 


0.4% 


M > 4.0 


70 


1 


1.00 


0.077 


59.5% 



